International Journal of Advanced Engineering Research and Science (IJAERS) 
httos://dx.doi.ora/l 0.22161/iiaers.5.7.26 


[Vol-5, Issue-7, July- 2018] 
ISSN: 2349-6495(P) / 2456-1908(0) 


Approach of economic-emission load dispatch 
using Ant Lion Optimizer 

Jorge de Almeida Brito Junior 1 , Manoel Henrique Reis Nascimento 2 , Carlos 
Alberto Oliveira de Freitas 3 , Jandecy Cabral Leite 4 , Tirso Lorenzo Reyes 

Carvajal 5 


Research department, Institute of Technology and Education Galileo of Amazon - ITEGAM, Manaus, AM, Brazil 

Email: jorgebritojr@gmail.com 

2 Research department, Institute of Technology and Education Galileo of Amazon - ITEGAM, Manaus, AM, Brazil 

Email: hreys@itegamorg.br 

3 Research department, Institute of Technology and Education Galileo of Amazon - ITEGAM, Manaus, AM, Brazil 

Email: carlos .freitas @itegamorg.br 

4 Research department, Institute of Technology and Education Galileo of Amazon - ITEGAM, Manaus, AM, Brazil 

Email: jandecy.cabral@itegamorg.br 

5 Research department, Institute of Technology and Education Galileo of Amazon - ITEGAM, Manaus, AM, Brazil 

Email: tirsolrca@gmail.com 


Abstract — To solve the problem of the economic 
emission load dispatch (EELD) is necessary minimize the 
total cost of fuel consumption and carbon emission. In 
this study is applied the ant lion optimizer (ALO) to this 
problem. The cost function and emission function with 
their respective restrictions are being using. To present 
the results this proposal is applied in IEEE 30 bus system 
that consists of six thermal units. The results for this case 
study with the application of ant lion with all generators 
on with demand being met, the total fuel cost is 
48915.36652 ($/h). The results this method can be 
compared with another metaheuristic algorithms and 
helps the plant operators in the decision making of 
preventive maintenance. 

Keywords—Ant lion Optimizer, EELD, Power Plants. 

I. INTRODUCTION 

The Thermal Power Plant (TPP) operation is dependent 
upon incineration of fossil fuel which generates sulfur 
dioxide (SO 2 ), carbon dioxide (CO 2 ) and nitrogen oxides 
(N0 X ) which create atmospheric pollution. Reduce the 
emission level and total cost of generation and at the same 
time accomplishing the demand for electricity from the 
power plant is is the goal of economic emission load 
dispatch (EELD). To solve the EELD problem is 
necessary minimize the total cost of fuel consumption and 
carbon emissions (De, Das, Mandal, & Mandal, 2018; 
Moraes, Bezerra, Moya Rodriguez, Nascimento, & Leite, 
2018). The problem is formulated as a multiobjective 
economic emission load dispatch (EELD) problem in 
which both the objectives (emission and economy) have 


to be minimized (Chopra, Kumar, & Mehta, 2016). This 
is a complex problem to solve because of its large size, a 
nonlinear objective function and a wide number of 
restrictions (Bhattacharya & Chattopadhyay, 2010). 
Various evolutionary, heuristic and meta-heuristics 
optimization algorithms have been developed such as: 
Grey Wolf Optimization (GWO) (Chopra et al., 2016; 
Hong, MH, & Mohd Rusllim, 2014), non-dominated 
sorting genetic algorithm (NSGA-II) (Basu, 2008; Moraes 
et al., 2018), hybrid genetic algorithm (Thenmozhi & 
Mary, 2004), Tabu Search Algorithm (Ii, Yang, Tseng, 
Wang, & lim, 2018), Simulated annealing (Junior, 
Nunes, Nascimento, Rodriguez, & Leite, 2017; Ziane, 
Benhamida, & Graa, 2017), Neural Networks (Deng, He, 
& Zeng, 2017), Harmony Search Algorithm (El Ela, El- 
Sehiemy, Shaheen, & Shalaby, 2017), particle swarm 
optimization (De et al., 2018), Differential Evolution 
(Jebaraj, Venkatesan, Soubache, & Rajan, 2017), Ant 
Colony Optimization (Zhou et al., 2017), Biogeography- 
Based Optimization (Ma, Yang, You, & Fei, 2017), 
genetic algorithm controlled by fuzzy logic (Song, Wang, 
Wang, & Johns, 1997). 

This research use the emission function and economic 
function in the multiobjective optimization ALO, with 
restrictions. 

II. MATERIAL AND METHODS 

To solve a problem of EELD, two important objectives in 
an electrical power system must be considered; they are: 
environmental, and economy inpacts (Basu, 2014). 

2.1 Economic Load Dispatch 
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The fuel cost is considered as an essential criterion for 
economics analysis in ELD. The most simplified cost 
function of each generator can be assumed to be 
approximated by a quadratic function of generator power 
output Pi (Ghosh, Chakraborty, Bhowmik, & 
Bhattacharya, 2017; Jebaraj et al., 2017): 

N 

FI (Pi) = (a, + b t P t + c^ 2 ) %/h (1) 

i =1 

where au hi and a are the fuel cost coefficients of the ith 
unit generating, N the number of generators and Pi the 
active power of each generator. Fig. 1 illustrates the fuel 
cost curve without valve-point effects and emissions. 



Fig.l: Fuel cost and emission function of the thermal 
generator. 

Source: (Gitizadeh & Ghavidel , 2014) 


2.2 Economic Emission Dispatch 

Emissions can be represented by a function, that links 
emissions with power generated by each unit. The 
emission function in kg/h, which normally represents the 
emission of S02 and NOx, is a function of the power 
output of the generator, and it can be express ed as follows 
(Swain, Sarkar, Meher, & Chanda, 2017): 

N 

F2(P l ) = Y J (d i + e l P i +/^, 2 ) kg/h (2) 

i-1 

Where du ei and fi are the emission coefficients of the ith 
unit generating, N the number of generators and P, the 
active power of each generator, from the TPP. 

2.3 Economical load dispatch constrains 
2.3.1 Equality power balance constraint 

The real power of each generator is limited by the lower 
and upper limits. The following equation is the equality 
restriction of power balance (Rizk-Allah, El-Sehiemy, & 
Wang, 2018). 

n 

£ P t - P D - P L (2) 

i= 1 


where Pi is the output power of each i generator, P D is the 
load demand and P L are transmission losses, in other 
words, the total power generation has to meet the total 
demand P D and the actual power losses in transmission 
lines P L (Dewangan, Jain, & Huddar, 2015). 

n 

Y J P i =P D -P L ( 3 ) 

i = l 

The calculation of power losses P L involves the solution 
of the load flow problem, which has equality constraints 
in the active and reactive power on each bar as follows 
(Nwulu & Xia, 2015): 


n 

p l =Yj BiP ? (4) 

i =1 

A simplification is applied to model the transmission 
losses, setting them as a function of the generator output 
through Kron’s loss coefficient derivatives of the Kron 
formula for losses (Huang et al., 2018). 


n n 


n 




( 5 ) 


i =17 = 1 i = 1 

where Bij, Boi and Boo are the energy loss coefficients in 
the transmission network and n is the number of 
generators. A reasonable accuracy can be obtained when 
the actual operating conditions are close to the base case, 
where the B coefficients were obtained (Gitizadeh & 
Ghavidel, 2014). 

2.3.2 Production Capacity Constraint 

The power capacity total generated from each generator is 
restricted by the lower limit and by the upper limit, so the 
constrain is (De et al., 2018): 

P < P- < P (6) 

where 77 is the output power of the i generator, P m in.u is 
the minimal power of the i generator and Pmax.i , the 
maximal power of the i generator. 

2.3.3 Fuel Delivery Constraint 

At each time interval, the amount of fuel supplied to all 
units must be less than or equal to the fuel supplied by the 
seller, i.e. the fuel delivered to each unit in each interval 
should be within its lower limit F m in,i and its upper limit 
Fmax,i so that (Qu et al., 2018): 

F min.i ^ F im ^ F max.v 1 E N > m E M, (7) 

where Fi, m is the fuel supplied to the engine i at the interval 
m , Fi, m in is the minimum amount of fuel supplied to i 
generator and F max ,i is the maximum amount of fuel 
supplied to i generator. 

2.3.4 Optimization problem 

The multi-objective optimization problem is defined as 
follow: 

Minimize (P) = [FI(P),F 2 (P)] (8) 

where Fi(P), F 2 (P) are the objective functions to be 
minimized over the set of permissible decision vector P. 

2.3.5 Incremental fuel cost method 
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The incremental fuel cost can be obtained from the 
following equation (Tiwari, Dave, & Dwivedi, 2017): 

IC t = {l-CLi-Pgi + bi) %/MWh (9) 

where ICi is the incremental fuel cost at are the values of 
the different points of the actual curve of the incremental 
cost and hi are the values of the points on the 
approximated curve (linear) of incremental cost. P g i is the 
total power generation. The curve of incremental fuel cost 
is show in the following Figure. 2: 



Fig.2: Incremental Cost Curve of Power Generator. 
Source:(Nascimento, Nunes, Rodriguez, Leite, & Junior, 


2016) 


2.4 Ant lion optimization 

The Ant lion Optimizer (ALO) is a algorithm inspired 
by nature (Miijalili, 2015). The ALO algorithm mimics 
interaction between antiions and ants in the trap. To model 
such interactions, ants are required to move over the 
search space, and antiions are allowed to hunt them and 
become fitter using traps. Since ants move stochastically 
in nature when searching for food, a random walk is 
chosen for modelling ants’ movement as follows [28]: 

X(t) = [0 ,cumsum (2r(t 1 ) — 1), cumsum(2r(t 2 ) 

— 1);..., cumsum(2r(t n ) — 1)] (10) 


where cumsum calculates the cumulative sum, n is the 
maximum number of iteration, t shows the step of random 
walk (iteration in this study), and r(t) is a stochastic 
function defined as follows (Trivedi, Jangir, & Parmar, 
2016): 




(u) 


(1 if rand > 0.5 
1 if otherwise 
where t shows the step of random walk (iteration in this 
study) and rand is a random number generated with 
uniform distribution in the interval of [0, 1]. 

To keep the random walk in the boundaries of the search 
space and prevent the ants from overshooting, the random 
walks should be normalized using the following equation 
(Yao & Wang, 2017): 

(xf-at)x(4~cf) 

( Pi~ai ) 


$=■ 


■+cf 


( 12 ) 


where c\ is the minimum of i-th variable at t-th iteration, 
d\ indicates the maximum of i-th variable at t-th iteration, 


a t is the minimum of random walk of i-th variable, and b t 
is the maximum of random walk in i-th variable. 

To simulate the trapping of ants the mathematical 
expression of the trapping of the ants to the ant lion’s pits 
is given by following equations (Trivedi et al., 2016): 
c l m — Ant — lion* — c t (13) 

d l m — Ant — lionf — d t (14) 

To construction of trap, the fittest ant lion is selected 
using the roulette wheel method. 

To simulate the sliding ants towards ant lions, the 
boundaries of random walks should be reduced adaptively 
as follows (Mirjalili, 2015): 
c t 

c t = j (15) 

f d ( 

d ( = j (16) 

where / = 10 w -, t is current iteration, S is the maximum 


number of iterations and w is a constant whose value is 
given by (Raju, Saikia, & Sinha, 2016): 
f 2 if t > 0.15 
3 if t > 0.55 

w = < 4 if t > 0.755 (17) 

5 if t> 0.95 
<6 if t > 0.955 

To catching the ants by ant lion and re-building the pit can 
be mathematically described as [28]: 

Antiion■ = Anti, if f(Antf) > f{Antlion-) (18) 

where Antlionj indicates the position of selected jth ant 


lion at tt h iteration and Antl shows the position of ith ant 
at tt h iteration, t shows the current iteration. 

Finally the last operator in ALO, that is elitism, calculated 
using roulette wheel as follows equation (Trivedi et al., 
2016): 


Antf ~ 



(19) 


where, R A = the random walk nearby the ant lion chose by 
means of the roulette wheel at t th iteration, R E = the 
random walk nearby the elite at t ih iteration, Ant\ = the 
location of / th ant at t ih iteration. 


2.5 ALO applied to EELD 

Initialize random walks on ants using Eq (10) and save 
generation scheduling of generating units as ant position 
using Eq (20) described below: 



Ant xl Ant 12 

... Ant lrd - 


^Ant ~ 

Ant 2 i Ant 2 2 ^^ 2,3 

... Ant 2d 

(20) 


Ant nl ... 

Ant nd . 

nxd 

where M Ant is the matrix for saving the position of each 


ant, Ant t j shows the value of the yth variable (dimension) 
of ith ant, n is the number of ants, and d is the number of 
variables. 
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For evaluating each ant (i.e., generating units), the 
following objective functions described in Eq. (1) and Eq 
(2) are utilized during optimization and following matrix 
stores the fitness value of all ants: 


f([Ant X i,Ant 12r - 


f([Ant 2i i,Ant 2 , 2l - 

■■Ant 2a ]') 

f(\Ant n i, Ant n 2 , • ■ 

•,Ant nra \) 


( 21 ) 


where M 0A is the matrix for saving the fitness of each ant, 
Anti j shows the value of jth dimension of ith ant, n is the 
number of ants, and/is the objective function. 

Save the optimal costand generation scheduling using 
Eqs. (22) and (23) described below: 


f(\AL llf AL 12f ■■■ j4£ad]) 

/([^2,li ^^2 r 2i ^2.d]) 

fi[AL n i t AL n 2, AL n ,$]) 


(23) 


where M 0AL is the matrix for saving the fitness of each ant 
lion, AL t j shows the jth dimension’s value of ith ant lion, 
n is the number of antiions, and/is the objective function. 

This solution comprises the number of generations of 
the system that will be optimized, which results in 
minimization of cost and emissions described in Eq (8) by 
fulfilling all constraints described in Eq (3), Eq (6) and Eq 
(7). 

Equation (8) are applied in the performance evaluation 
of the FEED until the optimum cost and emission is 
achieved. For inequality constraints, similar to any other 
techniques, when the solutions obtained for any iteration 
are out of boundaries, ALO chooses theboundaries values, 
while for equality constraint, when it is violated, the 
penalty factor of 1000 is implemented and embedded in 
the cost function as per Eq. (8). The algorithm will 
continue until the maximum iteration is met, and the 
optimum results are obtained. 

III. SIMULATION TESTS AND RESETS 
The power plant selected for the case study consists of six 
generating units with a load demand of 900 MW where 
generation limits, fuel cost and emission coefficients for 
case study is take from Ref (Lee & Darwish, 2008; 
Manteaw & Odero, 2012). 

The EELD problem simulated with the ALO algorithm, 
the systems of standard IEEE 30 bus systems have been 
taken into consideration (figure 3). 



Fig.3: Diagram of IEEE 30-bus test system. 
Source: (Lee & Darwish, 2008) 


The data of IEEE 30 bus test system to apply in ALO 
optimizer is presented in table 1, table 2 and table 3. 


Table. 1: Characteristic data of the generators from the 
case study power plants. 


Gen 

O 

($/MW 

2 h) 

bi 

($/M 

Wh) 

($/h) 

p . 

1 nun 

(MW 

) 

p 

1 max 

(MW 

) 


0.1524 

38.539 

756.79 

10 

125 

G1 

7 

73 

886 




0.1058 

46.159 

451.32 

10 

150 

G2 

7 

16 

513 




0.0280 

40.396 

1049.3 

40 

250 

G3 

3 

55 

2513 




0.0354 

38.305 

1243.5 

35 

210 

G4 

6 

53 

311 




0.0211 

36.327 

1658.5 

130 

325 

G5 

1 

82 

696 




0.0179 

38.270 

1376.2 

125 

315 

G6 

9 

41 

7041 




Source: (Manteaw & Odero, 2012) 


Table.2: Coefficients of emission of the 6 generating unit. 

Unit di e/ f 

G1 0.00419 0.32767 13.85932 

G2 0.00419 0.32767 13.85932 
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Source: Authors. 


G3 

0.00683 

-0.54551 

40.2669 

G4 

0.00683 

-0.54551 

40.2669 

G5 

0.00461 

-0.51116 

42.89553 

G6 

0.00461 

-0.51116 

42.89553 


Source: (Manteaw& Odero, 2012) 


Table.3: Loss coefficients (xlO 6 ). 


2022 

-286 

-534 

-565 

-454 

-103 

-286 

3243 

16 

-307 

-422 

-147 

-535 

16 

2085 

831 

23 

-270 

-565 

-307 

831 

1129 

113 

-295 

-454 

-422 

23 

113 

460 

-153 

-103 

-147 

-270 

-295 

-153 

898 


Source: (Manteaw& Odero, 2012) 


The results after running the simulation of the proposed 
ALO algorithm, are displayed in Tables 4.The simulation 
of the proposed ALO algorithm is tested in MATLAB 
R2016 to meet the demand of 900MW. 


The graphics with power, emission and cost are presented 
in figure 5, 6 and 7 respectively. 


Power generated 



1 2 3 4 5 6 

Generators 

Fig.5: Power generated in each generator. 
Source: Authors. 


Table.4: Coefficients of emission of the 6 generating 
unit. 


Ge 

ner 

ato 

r 

Power of 

each 

generator in 
Mw 

Emission of 
each engine 
in Kg/h 

Cost of each 
engine in $/h 

G1 

125 

8238.41182 

7956.60886 

G2 

92.7026704 

3382.26922 

5640.22656 

G3 

86.4365762 

4403.87432 

4750.48463 

G4 

151.819543 

23857.7694 

7876.38287 

G5 

240.571054 

64104.5546 

11619.7208 

G6 

229.179195 

55417.1514 

11071.9428 

Tot 

al 

925.709039 

159404.0308 

48915.36652 


Source: Authors 


The graphics with pareto front of costs versus emissions 
and the using all generators is presented in fig. 4. 


x10 s ALO OPTIMIZATION 



Fig.4: Pareto front of cost vs emission. 
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jET 6 
Si 
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— 
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0 

1 2 3 4 5 6 

Generators 

Fig.6: Emission in each generator. 
Source: Authors. 
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Fig. 7: Cost of each generator. 
Source: Authors. 
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iy! conclusion 


The ant lion optimizer is successfully applied to a 30 bus 
test system, to solve the EELD problem, so now it is 
possible to use these results to compare with other 
techniques that apply to this same IEEE bus test system 
system. This application can also help workers to operate 
more efficiently the generators in a power plant. 
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